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The Fano-Kondo effect in zero-bias conductance is investigated based on a theoretical 
model for the T-shaped quantum dot. The conductance as a function of the gate voltage 
is generally characterized by a Fano asymmetric parameter q. With varying temperature 
the conductance shows a crossover between the high and low temperature regions compared 
with the Kondo temperature Tk- two Fano asymmetric peaks at high temperatures and the 
Fano-Kondo plateau inside a Fano peak at low temperatures. Temperature dependence of 
conductance is calculated numerically by the Finite temperature density matrix renormal- 
ization group method (FT-DMRG). 
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When a localized state is coupled with a continuous spectrum phase relation of the in- 
terference is important and this phenomenon is known as the Fano effect. 1 In many cases, 
Coulomb interaction in the localized state cannot be neglected and it leads to the interesting 
many-body problem of the Kondo effect. However in usual Kondo effects of magnetic impu- 
rities the former aspect is trivial or not apparent and thus has not been discussed explicitly. 

Recently it has become possible to observe Kondo effects in transport properties through 
quantum dots (QDs). In the simplest geometry of a QD attached to source and drain leads, 
an embedded QD system, conductance reaches the unitarity limit 2e 2 /h at zero temperature 
for the single channel case in the Coulomb blockade region with odd number of electrons. 
Actually the unitarity limit due to the Kondo effect was observed experimentally. 2 

With the use of nano-scale technology it is possible to fabricate various geometries for 
QDs. The Fano effect combined with the Kondo effect for a single-level QD has been proposed 
theoretically in two geometries: a QD in a closed Aharonov-Bohm interferometer (AB-QD) 3,4 
and a side-coupled QD. 5 ~ 7 The latter geometry is a QD coupled with a quantum wire like a 
stub and we call it a T-shaped QD in the present paper. 

The line shape of zero-bias conductance as a function of the gate voltage is characterized 
by the Fano asymmetric parameter q. For the embedded QD, the conductance has a symmetric 
peak structure, which corresponds to infinite q. On the other hand conductance for the T- 
shaped QD sometimes shows a symmetric dip structure, which means q = 0. In this case 



1/9 



J. Phys. Soc. Jpn. Letter 

due to the Kondo effect the conductance becomes zero at T = and the system shows a 
perfect reflection rather than the resonance tunneling for the infinite q. This phenomenon is 
called as an anti-Kondo resonance in literature. 6 For an AB-QD, the Fano parameter takes 
a non-zero value. It can be controlled by the magnetic flux penetrating in the AB ring and 
can be a complex number. Experimentally, tunability of the Fano parameter in the AB-QD 
has been demonstrated successfully. 8 Fano effects with finite q are general phenomena in 
quantum transport and have been observed in various geometries not only in the AB-QD but 
also in single electron transistors 9 and T-shaped QDs. 10 However, so far, there has been no 
experimental report on indication of the Kondo effect for the cases with finite q except for the 
very recent experiment of the anti-Kondo resonance. 11 

Concerning theories for the T-shaped QD, the Fano-Kondo effect studied in ref. 5-7 has 
been limited to the q = case. In this letter we study the Fano-Kondo effect covering the whole 
range of q by considering an extended theoretical model for the T-shaped QD. In particular 
temperature dependence of the conductance is discussed for finite q and finite U. We note 
that the theories on the Fano-Kondo effect for the AB-QD have been limited to either the 
T = case 3 or the U = oo case. 4 It is also shown in this letter that the Fano-Kondo effect of 
the AB-QD can be described in the same way by the mapping to the extended model. 

We introduce the model for the T-shaped QD shown in Fig. 1 (a). Since we have introduced 
the energy level eo for the 0-th site the system may be considered as a T-shaped double dot 
system. 12 ' 13 The Hamiltonian of the model is 

U = -J2U ii+1 (c\ a Ci +lr7 + h.c.) 

in 

a 

-Vd^2(dlc 0a + h.c.) 

<j 

+e d d Ua + U d d\d^d\d b (1) 

a 

\ v ,when i = or -1, 
H,i+1 ■- \ {A) 
y t ,others., 

where each dot has a single level, eo and e^, and the on-site Coulomb interaction, Uq and 
Normalized conductance g = G j ' (jjf*) is given by 

xh (-?m d . , 3 , 



= J™ AoM(-Im[C?oM]) (-^£r) dw ' ( 4 ) 

where Ao(o>) = Ao^ 1 — {%) 2 with Ao = 2f 2 /i and Go(u) is the local Green's function of the 
0-th site. These formulae give conductance for the embedded QD when we cut the connection 
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to the d-dot by putting Vd = 0. In the following, we will concentrate on the T-shaped QD 
and Uo = is assumed unless otherwise stated. We also assume that the chemical potential 
is zero, which is the half-filling condition because the energy band of the bulk is — 2tcosfc. 

For the T-shaped QD model we can rewrite the conductance by Gd, because Gd and Go 
are related by 

CM = g (u)(l + v 2 d go(u)G d (u)). (5) 

The Green's function go is the local Green's function at the site for the case of Vd = and 
is given by 

00 M = (uj - e - 2vlg B (uj)) \ (6) 

where gs is the local Green's function at the i = ±1 site for the case of vq = 0, i.e. at the end 
point of the semi-infinite chain: 



g B (uj) = / dz — . 7 

J-2t U ~ Z + 10+ 

Since we can expand go(uj) as go{oj) c± 5o(0) + 0o(O) 2 (l + ^o/^) w + 0(uj 2 ), we can use go{oj) = 
5o(0) at low temperatures, T <C A . With using the limiting value go(uj) = go(0) = — (eo — 
iAo) -1 , the transmission probability at T <C Ao and t is written as 

1 2q 
1 + g2 ~ TTq' 2 

i 2 -I 



T{co) = - « A,Re[G d (a;)] 



1 + 



rAdlmfGdH], (8) 



where the Fano parameter q is 



and the width of the d-level is given by 



^An 



A ^ " M (10> 



^5 

It should be noted that in the present model the conductance has a contribution from 
the real part of Gd(co) for non-zero q. Especially at \q\ = 1, only Re[Gd(w)] contributes 
to the conductance. For other two limiting cases, q = and \q\ ~ oo, the conductance is 
determined solely by Im[G(j(a;)]. The former case, q = 0, is called anti-resonance case where 
the transmission probability is given by 

T(u>) = 1 - (-AdlmGdH) . (11) 

The latter case, \q\ ~ oo, is written as 

T(u) ~ -A d lmG d (oj). (12) 
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Fig. 1. (a) The model for the T-shaped quantum dot. (b) Mapping to two semi- infinite chains of the 
symmetric Sj and antisymmetric a, sites, Sj = (c, + c_ i )/v / 2 and a, = (q — c_i)/\/2. 



It is worth to mention that conductance in the embedded QD corresponds to the \q\ = oo 

case, 

T{lo) = -A ImG (cj). (13) 

Equivalence between the embedded QD and the \q\ = oo case of the T-shaped QD may be 
understood in the following way. In the limit of \q\ — > oo, we can write A^ = | ^^7^^ ; which 
shows the second order process plays the role of effective hopping between the d-dot and the 
leads. 

It is instructive to consider a mapping from the Hamiltonian of the T-shaped QD to two 
semi-infinite chains. The mapping is obtained by considering the symmetric Sj = (cj + c_j)/\/2 
and the antisymmetric ai = (q — c_j)/\/2 states, Fig.l (b). From the mapping it is readily 
seen that G,i{oj) with ^/2vq = t and eo = is identical to Gq(uj) with v,i = when U$ = Ud- 
However the conductances are different as shown in eqs.(ll) and (13). The Kondo resonance 
in the embedded QD corresponds to the anti-Kondo resonance in the T-shaped QD at eo = 0. 

It is easy to calculate the non-interacting case. Since the Green's function is calculated as 

GdM = ^~~2 — rr> (14) 

the conductance at T = is given by 

' = rh^- < 15 > 

1 + q z e z + 1 



with 

For interacting case at zero temperature we introduce the phase e of the Green's function 
at lo = as 

G d (0) = t^^-t, (17) 
A d -e + i 

where is defined by eq.(10). This is valid for any finite Ud, which is an essence of the 
Kondo effect. The conductance for general Ud at the zero temperature limit is given by the 
same formula, eq.(15). These formulae of g and q are the same as the non-interacting case. 
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What is different from the non-interacting case is e which contains all effects of interaction at 
zero temperature and defines the phase shift ^AiV: 

e = cotgAiv), (18) 

where AiV is the total number of up and down electrons induced by the impurity, i.e. < 
AiV < 2. The definition is AiV = (N) — (-/V)b u ik> where (iV)buik ^ s defined by the Hamiltonian 
without the d-site. The present relation is nothing but the Friedel sum rule. 

To study the conductance of the T-shaped QD numerically, we use the Finite Temperature 
Density Matrix Renormalization Group (FT-DMRG) method. This method has an advantage 
that we can obtain results for different impurities in the same host one-dimensional system only 
by taking the expectation value concerning the impurity part. Since a thermal Green's function 
can be calculated with the FT-DMRG, the local Green's function Gd(u>) can be obtained from 
the thermal Green's function G^(r) with numerical analytic continuation, where we use both 
the Pade approximation and the Maximum Entropy Method (MEM). For the FT-DMRG, the 
trotter number M < 200 and the residual bases m < 145 are used. For further details of this 
method we refer the readers to ref. 14. 

At high temperatures conductance obtained by the FT-DMRG shows typically two struc- 
tures at V g = ±U/2, i.e. two dip structures for q = (Fig. 2) and two Fano line shapes for q / 
(Fig. 3). In these figures, we plot also conductance obtained by the Zubarev approximation: 

G d (u) = (FW-vjgoiu))- 1 (19) 

FM = (l^ + -^L-)~\ (20) 
V u;-e d u-e d -UJ 

where (n d ) = (n^) = (ridj.) is the number of local electrons. In the figures (n d ) obtained 

by the FT-DMRG is used in order to avoid further approximation to calculate (rid). These 

figures show perfect agreement at U d = data between the FT-DMRG and the Zubarev's 

approximation which is exact in the noninteracting case. The Zubarev's approximation shows 

reasonable agreement with the FT-DMRG results at larger Ud- Characteristic features of 

conductance at high temperatures may be understood from the Zubarev's approximation. 

However, it should be noted that the Zubarev's approximation is not a controlled scheme 

even at high temperatures, which is the origin of the differences between the FT-DMRG 

results and the Zubarev's approximation observed in Figs. 2 and 3. 

Both the Zubarev's approximation, eq.(20), and the zero temperature expression, eq.(18), 

are written by the number of electrons, (n d ) and AiV. Let us suppose that temperature 

dependence of the number of electrons can be ignored and (n d ) is given by Fig. 4(a). At 

high temperatures we see two Fano structures at around — U/2 and U/2 on the background 

conductance of j^p, Fig. 4(b). At low temperatures we see development of a new Fano-Kondo 

2 

plateau of j^p, Fig. 4(c). At the same time the structure at —U/2 (U/2) becomes a simple 
peak (dip). 
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Fig. 2. Conductance of T-shaped QD for Ud — 0, 1, 2, 4 as a function of gate voltage V a = + U/2 
with q = 0, A /t = 1 and A d /t = 0.25 at T = 0.15. 
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Fig. 3. Conductance with q = —0.8 and Ad/t = 0.15. Other conditions are the same as Fig. 2 



Numerical results of temperature dependence of conductance are plotted in Fig. 5 for q = 
and in Fig.6 for q = —0.8. Differences between the results obtained by the Pade approximation 
and the MEM may be considered as typical errors of the numerical calculations. In Fig. 6, at 
lower temperatures only the results by the MEM are shown due to the numerical instability of 
the Pade approximation. These figures show smooth crossover from high temperature regime 
to the zero temperature limit. In these examples relatively large is used in order to achieve 
the low temperature regime within the present numerical scheme. Actually, the q = anti- 
Kondo resonance was recently observed by Katsumoto et al. 11 in a system with relatively 
large A^/U, the results of which are qualitatively consistent with Fig.5. 

At zero temperature, conductance is given by AN due to the Friedel sum rule, eq.(18). 
AiV is obtained from the numerical differentiation of the Free energy by changing the chemical 
potential /j,, (N) = —5F/5fj,. Conductance obtained from AN evaluated at T = 0.02t is shown 
in Figs. 5 and 6. In Fig.6 the conductance obtained by A7V extrapolated to T = is also 
shown. The difference is so small and cannot be distinguished. The Fano-Kondo plateau at 
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Fig. 4. (a) Number of electrons as test function for (rid) = AN/2, (b) The conductance of the 
Zubarev's approximation with q = —0.8 and (rid) defined in (a), which shows characteristic be- 
haviours in the high temperature regime, (c) The conductance in the low temperature regime 
given by the sum rule from AN/ 2 defined in (a). 
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Fig. 5. Numerical results of conductance as a function of V g for q = 0, Ud = 4 and Ao/t = Ad/t = 1 
at various temperatures calculated by the FT-DMRG with the MEM and the Pade approximation. 
go is the conductance of the background defined at \V g \ = oo limit. 



is not clear for this set of parameters. The main reason is relatively large A^/Ud- Another 
reason is that for a finite q the contribution from Re[Gy has a finite slope around V g = 0. This 
effect becomes important at q ~ ±1. We note that the center of this plateau where AiV = 1 



is shifted from V g = 0. 



Finally we comment on the AB-QD. The simplest tight binding model for this case is 
shown in Fig. 7 (a). For T <t the transmission probability is given by eq.(8) with 



1 



w 



2w 



(21) 



The model for the AB-QD may be mapped to two semi-infinite chains by introducing the 
symmetric, Sj, and the antisymmetric, ctj, sites. From the mapping it is readily seen that 
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Fig. 6. Numerical results of conductance for q = —0.8, Ud = 4, A /t = 1 and Ad/t = 0.61. At low 
temperatures, T — 0.06 and 0.1, only the results by the MEM are plotted due to the numerical 
instability of the Pade approximation. Other conditions are the same as Fig. 5. go is the conductance 
of the background defined at \V g \ = oo limit. 




42 Vd 

Fig. 7. (a) The model for the Aharonov-Bohm quantum dot. (b) Mapping to two semi- infinite chains 
by the symmetric Si and antisymmetric a,i states, Si — (ci + c_i)/V2 and <n = (ci — c-i)/^/2. 



Gd{oj) in the present model is essentially the same as the T-shaped QD with the identification 



Ad 2' ^ Z1 > 



where Ao = t in the present model. Therefore results discussed in the present paper are valid 
also for the AB-QD. 
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